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i- 1 Abstract 

Monopoles in topologically massive gauge theories in 2+1 dimensions 
with a Chern-Simon mass term have been studied by Pisarski some years 
r*| ago. He investigated the SU(2) Yang-Mills-Higgs model with an additional 

Chern-Simon mass term in the action. Pisarski argued that there is a 
£S| monopole solution that is regular everywhere, but found that it does not 

possess finite action. There were no exact or numerical solutions being pre- 
sented by Pisarski. Hence it is our purpose to further investigate this solu- 
fN| tion in more detail. We obtained numerical regular solutions that smoothly 

interpolates between the behavior at small and large distances for different 
values of Chern-Simon term strength and for several fixed values of Higgs 
^ field strength. 



1 Introduction 

X 

The study of 2+1 dimemsional gauge theories with a Chern-Simons (CS) term 
added to the Lagrangian density was first proposed by Deser et al. PQ. They 
were motivated by the connection of 2+1 dimensional gauge theories to the high 
temperature limit of four dimensional models. They observed that the massless 
gauge field which is spinless becomes massive with spin 1 and the quantum theory 
is also infrared and ultraviolet finite in perturbation theory when interacting with 
fermions. Furthermore topology also leads to a quantization condition on the 
coupling constant-mass ratio, ^jf = n, where the integer n = 0, ±1, ±2, . . . , in the 
non-Abelian gauge theories. 



*to be submitted for publication 
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Various non-Abelian gauge theories were studied by various authors in the 
1980's. One of these is the work by Pisarski [2J. In his work on the 2+1 di- 
mensional Yang-Mills-Higgs (YMH) theory with the CS term in the Euclidean 
three space, he noticed that the one-monopole solution analogous to that of the 
't Hooft-Polyakov one-monopole solution of the SU(2) Yang-Mills-Higgs (YMH) 
field theory |3] exist. This famous 't Hooft-Polyakov one-monopole solution with 
nonzero Higgs mass and self-interaction is a numerically spherically symmetric, fi- 
nite energy solution in 3+1 dimensions. Contrary to Polyakov [4J, who has shown 
that in 2+1 dimensions without the CS term, a dilute plasma of monopoles pro- 
duces confinement of external U(l) charges, Pisarski found that it is the monopoles 
that are confined and not the external electric charges with the CS included in the 
theory. He also remark that the analogous monopole solution that he had found, 
has no counterpart in 3+1 dimensions. Pisarski argued that there is a monopole 
solution that is regular everywhere, but found that it does not possess finite action. 
There were no exact or numerical solutions being presented by Pisarski. Hence 
it is our purpose to further investigate this solution in more detail. We obtained 
numerical regular solutions that smoothly interpolates between the behavior at 
small and large distances, for different values of Chern-Simon term strength and 
for several fixed values of Higgs field strength. 

In the next section, we will discuss briefly the 2+1 dimensional YMH theory 
with the CS term. In section 3, we discuss the magnetic one-monopole ansatz and 
in section 4, we present our numerical one-monopole solutions. We end with some 
remarks in the last section. 



2 The SU(2) Yang-Mills-Higgs Theory with the 
Chern-Simons term 

The SU(2) YMH Lagrangian with non vanishing Higgs potential is given by 

11 1 ^ 

Cymh = —F%F a ^ + -D^ a D^ a - -A($ a $ a - ^-f . (1) 

Here the Higgs field mass is /i and the strength of the Higgs potential is A. The 
vacuum expectation value of the Higgs field is £ = fi/VX. The YMH Lagrangian 
([I]) is gauge invariant under the set of independent local SU(2) transformations at 
each space-time point. The covariant derivative of the Higgs field and the gauge 
field strength tensor are given respectively by 

D^ a = d^ a + ge abc A b ^ c , 
F; v = d,Al-d u Al + ge abc AlA c u . (2) 

In the non-Abelian SU(2) YMH theory, the inclusion of the CS term leads to 
the Lagrangian density which is given by 

£ = £-ymh + £-cs (3) 
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where the YMH Lagrangian density is given as in Eq.([2j) and the CS Lagrangiann 
density is given by 

£cs = \t<r* (Aid u A a a + \a%M^ . (4) 

The CS parameter £ and the square of the gauge couping constant g have the 
dimension of mass. The CS parameter £ is real in Minkowski space and is replaced 
by — i£ in the Euclidean space. The SU(2) group indices a, b, c run from 1 to 3 and 
the space-time indices [i, v, a — 1, 2 and in the Minkowski space and 1, 2 and 3 
in the Euclidean space. The Minkowski space metric used here is g^ v = (+ H — ). 

The equations of motion that result from the above Lagrangian density Eq. ^ 
are gauge covariant, 

D^F; U + ^e vm F^ a = e abc $ b D u <S> c , 
where D»F^ = + g^A^F^, (5) 

and D^D^ a = -A$ a ($ 6 $ 6 - ^-), (6) 

A 



even though the Lagrangian density (|3j) is not gauge invariant. Under a finite 
gauge transformation uj(x), where the gauge potentials are changed by, 

A„ ->■ A^ = ujA^uj- 1 - (d^u- 1 , (7) 

the action S = J d 3 xC transforms as 

S -)■ 5" = [d 3 x e^TtdJA^d^) 

g 2 J 

+ f d 3 x e^Triu^daiouj^dpuu^d^u). (8) 
6g J 

In order for the first surface integral term in Eq.Q to vanish, u(x) must tend to 
the identity at time and spatial infinity, that is u(x) — > I. By explicitly writing 
oj(x) to be \I\ 

u(x) = exp{-6 a (x)T a ) , (9) 

g 2 



W{u) = —— - I d 3 x e Q/?7 Tr (u^d^u^dauu^d-u) 

24-7T 2 J ^ ' 

= ^Jd 5 x e^e^d a Ud,9%9^ 2 (l + Sk f\ } , (10) 



where # 2 = # a a and W(u) which is an integer is the winding number of the gauge 
transformation u. Here T a = ^ is the SU(2) gauge transformation generators and 
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a a are the Pauli matrices. The winding number W(u) characterizes the homotopic 
equivalent class of u(x) and vanishes only when u(x) can be continuously deformed 
to I. The quantization of W(oj) is valid both in the Euclidean as well as in the 
Minkowski space because it is a metric-independent coordinate invariant. 

Hence the action S is not gauge invariant but the exponential of the action, 
exp (i J d 3 xC) has to be gauge invariant. Therefore 87r 2 £/g 2 must be equal to 2nn 
giving 

4<// = n, n = 0,±l,±2... (11) 

The same is true in the Euclidean space as the requirement that exp / d?xC be 
gauge invariant still holds as £ — > — i£ in the Euclidean space. 
The dual field which is defined by 

* F »a = e ^i^ (12) 

satisfies the Bianchi identity D^*F® = 0. The usual physical quantities which 
characterise the solutions are given as follows: 

(13) 

(14) 
(15) 

(16) 
(17) 
(18) 

where i,j = 1,2 and r] a (x) is a unit vector in the internal group space that trans- 
form gauge covariantly. In all our discussion we set the gauge coupling constant g 
to unity. 

3 The Magnetic One-Monopole Ansatz 

The magnetic one-monopole ansatz used in Euclidean space is [2] 

gAl = ^(l-r 1 (r))(^-0 a ^) + ^r 2 (r)(O M + a ^) + f%A(r), 

g <S> a = $( r )f a . (19) 



Magnetic field: 




1 

= "2 e 




Electric field: 


E? 


= K 




Energy-momentum tensor: 




= F^ a F» a + g^L YM 


Angular momentum: 


J 






Magnetic flux: 






d 2 xB b r] b 


Electric charge: 


Q a 




d 2 x{d i E\)rf > 



where r = J x\ + x\ + x\ and the spatial spherical coordinate orthonormal unit 
vectors are defined by 

ffi = sin 9 cos0 5 M i + sin 9 sin0 <5^ 2 + cos 9 5^, 
9^ = cos 9 cos0 8^1 + cos 9 sin0 5^2 — sin 9 5^, 
M = -sin0 + cos0 5^2, (20) 
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The ansatz ( 19 ) is the most general ansatz that is invariant under combined rota- 

(21) 



tions of isospin and space-time and under the gauge transformation 

u{x) = exp (-0° (x)T a ) , 6 a = f{r)f a , 
the gauge potentials transformed as in Eq.(j7]) giving 



A 



n cos f + r 2 sin /, 

-Ti sin / + r 2 cos /, 
df 



A 



dr 



(22) 



The total action S in the three dimensional Euclidean space is given by 

S : 

Symh 



Symh + Scs 
1 



d x 



pa paixv + \ D H$a D $ a _ I 

4 2 4 



A^\2 

A ' 



roc ( 1 

I dr{{T[ + AT 2 y + {T' 2 -AT i y + —{l-Tl 



+ 



47T 

47r r 00 



r. 



2\2 



5, 



CS 



g 2 Jo 



d 3 x 



dr < (tI + tI) 



A 2 $ 4 

2 4 



1 



Air 



g 2 Jo 



POO . - 

y rfr^{r(r 2 + r^(l-r 1 )-A(l-r 1 2 -r 2 2 )}. 



Under the gauge transformation of Eq.(22), the action changes by 



Scs = S C s + {r 2 (cos / - 1) - n sin / + /} 



(23) 



(24) 



The gauge transformations for which /(0) = 27m i and /(oo) = 27m 2 , where ni 
and n2 are integers, do not change T\ and t 2 at r = and r = oo. However the 



action (23) changes by i^^-n, n=integer. Since exp(zS') is invariant, we have 



the quantization condition -^£=integer, that is Eq.(ll). Hence the presence of 
monopoles do not change the quantization condition. 



Substituting the ansatz (19) into the equations of motion ([5| and ^ and 
letting T\ — r, and r 2 = 0, the equations of motion (|5| and ^ reduced to only 
two coupled differential equations 



,-r(r 2 



D- 



e /r 4 -i' 



$ 2 r = 0, 



r" 



with A = — — . 



(25) 
(26) 
(27) 



Pisarski [2] claimed that there is a regular U(l) monopole solution of the above 



equations of motion (25) - (27). However he did not solved these equations of 
motion. In this paper, we solved these equations of motion numerically and prove 
that these monopole solutions do exist. 

These regular solutions are monopole solutions as the nonzero Higgs field points 
in the r a direction in internal space. This is in accordance to the electromagnetic 
field as proposed by 't Hooft [3], 

7 — <f>" 



P = |> a F° - e abc § a DA h D v & 



= d^A u -d u A^-e abc &d^ b d u &, (28) 
where A^ = ^ a A^, the Higgs unit vector, $ a = $ a /|$|, and the Higgs field mag- 



nitude |$| = \J $ a $ a . Hence the electromagnetic field here according to Eq.(28) 
is 

i^=(^-*A)J ( 29 ) 

which is the magnetic field of a point monopole. 

The boundary conditions that we consider for our numerical calculations are 

t ->■ 1, A -> 0, $ -> r, as r^O, (30) 
i2 -\ i/4 



$->C = y^- as r -»• oo. (31) 



Pisarski cannot find an exact solution with the above boundary conditions (30) 



and (31). However, he did mentioned that surely there is a regular solution that 
smoothly interpolates between the behavior found at large and small r [2] . He also 
mentioned that while this solution is regular, it does not possess finite action as 
the total action, 

4tt fR( /2 (l-r 2 ) 2 /2 e\ r 2 . 2 _ 2 /i 2 2 A 2 4 1 
S = — I { t' 2 + ; — — I — + \ ) + — $ /2 + $ 2 r 2 - ^r 2 ^ 2 + -r 2 $ 4 |> rfr 

! ' '^T^T- K (32) 



c/ 2 7o [ 4 \r 2 r 2 y ' 2 " 2 ' 4 

# 2 Jo 4 



is proportional to aR where is the radius of space-time and a is the string tension 
for a monopole-antimonopole pair and is given by 

(33) 

The plot in Figure 1 shows how the string tension a varies with the CS parameter 
£ in units of when £ = 1. Hence for small CS parameter £ <C £, cr and 
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Figure 1: The plot of the string tension a versus the CS parameter £ in units of 
assuming that £ can be plotted continuously. 



it increases proportionally with £. However for large CS parameter £ ^> £, cr — >• 
^fC 2 - Therefore the presence of the CS term not only causes all gauge fields, 
including the photons, to be massive but also causes monopole and antimonopole 
to be tightly bound together. However as argued by d'Hoker and Vinet [5], since 
there is no long range correlations, topologically massive theories do not confine 
external charges but instead monopoles are confined. 



4 The Numerical One-Monopole Solutions 



In this paper the numerical calculations were performed using the Maple 12 and 
MatLab 2007 softwares. The two second order equations of motion d5| - (|6]) are 
reduced to two differential equations (25) - (26) with the ansatz (19). These 



two differential equations are then transformed into a system of nonlinear equa- 
tions by using the finite difference approximation. They are then discretized into 
50 divisions in the x coordinate where x = -4r is the so-called finite interval 
compactified coordinate. By considering the x coordinate, the derivative with re- 
spect to the radial coordinate is then replaced accordingly by 4- — > (1 — x) 2 j|, 
A — >■ (1 —x) A -^2 — 2(1 — x) 3 -^. The system of nonlinear equations are then solved 
numerically with the Gauss-Newton method by providing the solver with good 
initial guess. 

When the expectation value of the Higgs field is ( = 1, and A = 0, |, 1, and 10, 



we solved the two differential ( 25 ) - ( 26 ) numerically for r and $, for different values 
of CS parameter £ = 0,0.1,0.5,1 and 10. These numerical solutions are shown 
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Figure 2: (a) The plot of $(r) versus r and (b) the plot of r(r) versus r for different 
values of CS parameter £ = 0, 0.1, 0.5, 1 and 10, when A = 0. 



graphically in Figure ([2]) - ([5]). We can conclude from the graphs of the numerical 
solutions that with increasing values of the strength of the Higgs potential A, the 
Higgs field approaches its expectation value £ = 1 faster. However the effect of the 
CS parameter slows down the rate of $ approaching (. Similarly with increasing 
values of the strength of the Higgs potential A, the function r — > faster. The 
asymptotic value of r at large r, that is r^, however now varies with the CS 

parameter £ as r M = { f2 £ * 2 } / . 



5 Remarks 



As mentioned by Pisarski [2] there are four solutions to Eq.(25) at large r, as 
r — > ifoo, iiroo. However the imaginary solutions can give negative values of the 
string tension a and hence are not physical. Also purely imaginary gauge potentials 
A® will also give negative action density and may not represent a physical vacuum 
state. 

We have found the regular numerical one-monopole solutions of the massive 
topological gauge theory whose existence was predicted by Pisarski [2] and he in- 
terpreted the string tension a as the tension between a monopole-antimonopole 
pair. However until now there were no monopole-antimonopole pair (MAP) solu- 
tions of the 2+1 dimensional YMH gauge field theory with the CS term. Hence 
our next task is to search for the MAP solutions in this theory. Similar to the 
MAP solution in 3+1 dimensional YMH gauge field theory, these MAP solution 
cannot posses spherical symmetry but only axial symmetry. 
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(b) 



Figure 3: (a) The plot of $(r) versus r and (b) the plot of r(r) versus r for different 
values of CS parameter £ = 0, 0.1, 0.5, 1 and 10, when A = 0.5. 

We would also like to note that in the 2+1 dimensional YMH field theory with 
the CS term, the one-monopole solution presented here is not a solition but instead 
it is an instanton and may represent some tunnelling events. 
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(b) 



Figure 4: (a) The plot of $(r) versus r and (b) the plot of r(r) versus r for different 
values of CS parameter £ = 0, 0.1, 0.5, 1 and 10, when A = 1. 




Figure 5: (a) The plot of $(r) versus r and (b) the plot of r(r) versus r for different 
values of CS parameter £ = 0, 0.1, 0.5, 1 and 10, when A = 10. 
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